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Abstract 

We present results of numerical studies of the Landau-Ginzburg 
dynamics of the order parameter in one-dimensional models inspired 
by the condensed matter analogues of cosmological phase transitions. 
The main goal of our work is to show that, as proposed by one of 
us 0, the density of the frozen-out topological defects is set by the 
competition between the quench rate — the rate at which the phase 
transition is taking place — and the relaxation rate of the order param- 
eter. In other words, the characteristic domain size, which determines 
the typical separation of topological defects in the new broken sym- 
metry phase, is of the order of the correlation length at the instant at 
which the relaxation timescale of the order parameter equals the time 
remaining to the phase transition. In estimating the size of topologi- 
cal domains, this scenario shares with the original Kibble mechanism 
the idea that topological defects will form along the boundaries of 
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independently selected regions of the new broken symmetry vacuum. 
However, it derives the size of such domains from non-equilibrium as- 
pects of the transition (quench rate), as opposed to Kibble's original 
proposal in which their size was estimated from the Ginzburg tem- 
perature above which thermally activated symmetry restoration can 
occur. 

PACS: 05.70.Fh, ll.15.Ex, 67.40.Vs, 74.60-w 

1 Introduction 

The expansion of the Universe inevitably leads to a drop of the temperature 
of the primordial fireball. This cooling provides natural conditions that are 
believed to precipitate a series of phase transitions in which "false" symmet- 
ric, high temperature phases are transformed into low temperature broken 
symmetry "true" vacuum. As the Universe undergoes such phase transitions, 
the selection of the low-temperature, broken symmetry phase can only occur 
locally, within the causally correlated regions. Kibble [0, [3| first noted that 
this symmetry breaking process may leave relics of the high energy phase 
which will be trapped by the topologically stable configurations of the bro- 
ken symmetry phase. In principle, such topological defects could be massive 
enough to leave an observational imprint in the cosmic microwave background 
and perhaps also produce the seeds needed for matter structure formation 

0' I" 

There are three principal types of topological defects |J. They differ by 
their dimensionality. Monopoles are point-like objects that would be disas- 
trous from the cosmological point of view. They have a tendency to dominate 
the matter content of the Universe. Inflation was invented in part to dilute 
them and thus to prevent observational inconsistencies. Membranes or do- 
main walls are two-dimensional, and they almost certainly could also lead 
to the same disaster as monopoles. Specifically, they would induce unac- 
ceptably large distortions of the cosmic microwave background. By contrast, 
one-dimensional topological defects, cosmic strings, could provide the source 
of density perturbations needed to seed structure formation ||. 

Symmetry breaking phase transitions which occur in condensed matter 
physics are described by theories which are formally analogous to those in- 
volved in the cosmological context (for a review see ||). Condensed matter 
phase transitions have, however, one crucial advantage: they can be studied 
in the laboratory. With this in mind, more than a decade ago one of us 
suggested [1], [7] that the cosmological mechanism for defect formation can 
be studied experimentally in the condensed matter context. A number of 
experiments have been carried out following that proposal. For instance, liq- 
uid crystal experiments demonstrated that copious production of topological 
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defects does indeed happen in symmetry breaking phase transitions [§, 0]. 
More recently, superfluid He 4 and He 3 experiments have also tested the topo- 
logical defect formation picture [lCj, O, 12|. Liquid crystal phase transitions, 
however, are of the first order type. In first order phase transitions, the pro- 
cess which is responsible for the appearance of the new phase is nucleation. 
Namely, small regions of the medium undergo thermal activation which takes 
them over the potential barrier separating "false" and "true" vacua. As a 
result, bubbles of a certain (critical) size appear and form seeds of the new 
phase. Eventually, through growth and coalescence of these bubbles, the 
new phase replaces the old one. A much more interesting and less under- 
stood dynamics takes place in second order (Landau-Ginzburg like) phase 
transitions. For this reason, He 4 and He 3 systems, for which the transition 
to a superfluid phase proceeds without nucleation, have recently attracted 
significant interest as examples of cosmological phase transition analogues. 

The defect formation scenario is based on two assumptions. The first 
assumption is of the qualitative nature: It asserts that regions of the broken 
symmetry phase, which are causally disconnected, must select the new low 
temperature phase independently. As a result, when a symmetry-breaking 
phase transition with a non-trivial vacuum manifold occurs simultaneously 
in a sufficiently large volume, topological defects will be formed. The second 
assumption is of the quantitative nature: It involves specifying the process 
responsible for the causal propagation of "signals" that allow the choice of 
the new vacuum to occur in a coordinated (rather than independent) fashion. 
This is the assumption that yields a prediction of the density of topological 
defects. Both assumptions are of course necessary 0, but while the first 
one is straightforward, the second one requires much more specific physical 
input. 

The original discussion of the scenario for defects formation was based 
on an idea similar to that of nucleation in first order phase transitions [^]. 
It was thought that the initial density of the topological defects is set at 
the so-called Ginzburg temperature T G , namely the temperature below the 
critical temperature Tq at which thermally activated transitions between the 
correlation-sized volumes of the new broken symmetry phase cease to occur. 
Below Tq-, the (free) energy barrier, separating "false" and "true" vacua, 
becomes prohibitively large for correlation-length sized thermal fluctuations. 
Under this approach then, the density of defects is set by the correlation 
length at T G . 

One of the key predictions of the original papers |], [/] on "cosmological" 
phase transitions in superfluids was that this thermally activated process 
does not decide the density of defects. Instead it was proposed that the 
characteristic correlations length is set by the dynamics of the order param- 
eter in the vicinity of the critical temperature Tq- This "non-equilibrium" 
scenario Jl|, [7j represents a significant change of the point of view, and yields 
a prediction for the density of defects which is rather different from the 
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"equilibrium" estimate based on the Ginzburg temperature. 



2 Symmetry Breaking 

In the vicinity of the critical temperature Tq of a second-order phase transi- 
tion, the contribution from the potential V to the free energy of the system 
is given by 

V= ~{\*\*-2o 3 e\*\ 2 + <T*€ i ), (1) 
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with A and a positive parameters and e = (T c — T)/T c the dimension- 
less temperature parameter. The last term in (ffl) has been added, so V = 
at the minimum of the potential for T < Tq- Here $ represents the 
order parameter, an abstract measure of the degree to which the symmetry 
in question has been broken. 

Above the critical temperature, the potential (II]) has a single minimum 
at = 0. As the temperature drops below T c , the shape of the potential 
changes; e becomes non-negative. Instead of a single minimum at |^| = 0, 
the potential has a minimum value V = at \^/\ = 0\ft. At the same time, 
the potential V develops a barrier centered at = with magnitude 

AV = V(0)-V(a v ^) = -Xcr 4 e 2 , (2) 
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"separating" the minimum values of the potential. 

Depending on the symmetry that is broken, the manifold of equivalent 
vacuum states for the order parameter will have different topologies, and 
thus describe a different set of topological defects [fL3| . For instance, in 



the case of domain walls, the order parameter is a real scalar field with 
discrete \1/ — > — ^ symmetry. The vacuum manifold in this case consists 
of the two values ^ = icry^e. Simple models of cosmic strings are, on the 
other hand, represented by a a complex scalar field with U(l) symmetry of 
phase transformations, \I/ — > e ia ^> . The degenerate vacuum states are given 
by \l/ = cr-y/ee* 61 , with 9 an arbitrary phase. Thus, the vacuum manifold for 
cosmic strings is the circle |\&| = cr-^/e in the complex ^-plane. 

As mentioned before, since all the points in the vacuum manifold are com- 
pletely equivalent, the choice of the expectation value of the order parameter 
is determined from the initial thermal (random) fluctuations through its dy- 
namics. A measure of the extent, in space and time, to which perturbations 
of the order parameter return to their equilibrium value is given by the corre- 
lation length £ and dynamical relaxation time r. Specifically, the correlation 
function of perturbations away from equilibrium behaves as 

{6*(x,t),6*{x + A,t)) ~exp(-|A|/0 (3) 
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and 

(8^(x, t), 5^(x, t + A)) ~exp(-|A|/r). (4) 

Furthermore, in the vicinity of the critical temperature Tc, the equilibrium 
values of £ and r are given by 

e = G>|e|-" (5) 

and 

r = To|e|-", (6) 

respectively, with /i and z/ > 0. That is, at the critical temperature (e = 0), 
the order parameter has a universal behavior characterized by the simulta- 
neous divergence of both £ and r. 



3 Freeze-out of Topological Defects 

Before presenting the results of the density of topological defects from our 
numerical study, let us review the original estimate of the density of defects 
put forward in This estimate was based on a linear quench, 

J-C Tq 

with t the time before (t < 0) and after (t > 0) the transition. (The idea 
behind the argument does not rely on this particular time dependence of the 
temperature parameter e, but it seems to be a likely approximation when 
narrow range of temperatures near critical Tc is involved, as it will be the 
case here.) In (0), tq represents the quench timescale. The other relevant 
timescale — which we shall call symmetry breaking timescale — is defined 
by |e/e| = \t\. It is then possible to identify during a quench three regimes 
by comparing the magnitude of the symmetry breaking timescale \t\ with the 
dynamical relaxation time r. 

For sufficiently large values of the temperature above T c (see Fig. |l|), the 
symmetry breaking timescale \t\ is larger than the dynamical relaxation time 
t. The system is in a stage where the order parameter is able to adjust to the 
changes implied by the quench in an essentially adiabatic fashion. In partic- 
ular, the non-equilibrium correlation length closely matches the equilibrium 
behavior given by Eq. (^|). 

Because of the critical slowing down implied by Eq. (§), at some time —t 
before the transition, the systems reaches the point at which the symmetry 
breaking timescale \t\ and the dynamical relaxation time r are comparable; 
namely 

= (8) 
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The value of the temperature parameter at this point is 

gM-*) = • (9) 

The corresponding value of the correlation length can be directly found by 
substituting Eq. (g) into Eq. (]5|), 

e(-i) = e P)^- (io) 

Within the [—t, t\ interval, the system is incapable of keeping up with the 
changes induced by the quench. In particular, the non-equilibrium correlation 
length cannot longer grow at the equilibrium rate implied by Eq. (|5]). Crudely 
speaking, the non-equilibrium correlation length "freezes out" at the value 
£(— t) in Eq. ([U]) until a time t, after the transition, when once again the 
dynamical relaxation time r becomes smaller that the symmetry breaking 
timescale \t\. That is, approximately, 

a-i) « e(t) = i (ii) 

The actual behavior of £(£) is probably more complicated (i.e. correlation 

will continue to grow even within the [—t,t\ interval), but £ does capture 
the size of the correlated region of vacuum at the time when its growth falls 
out of step with the equilibrium scaling. The bottom line is that during the 
interval [— t,t\, the order parameter, and therefore the correlation length, is 
to sluggish to keep up with Eq. (|5|). 

The proposal in [[[], [/J states that the frozen out correlation length £ 

yields the density of topological defects. That is, £ characterizes the size of 
independently selected domains of the new vacuum, so it determines the sep- 
aration (and therefore density) of topological defects. As mentioned before, 
the difference with the original Kibble mechanism j3| is the role that non- 
equilibrium dynamics play in determining the value of £. In that picture, £ 
was obtained entirely from the temperature Tq (Ginzburg temperature) at 
which fluctuations between the degenerate minima of the potential ceased to 
be possible; that is, the temperature at which 

k B T « | 3 AV (12) 

with AV the potential "barrier" in Eq. (|2|) separating degenerate minima. 

The reason why Tq does not appear to play a decisive role in the creation 
of topological defects is likely to be associated with the spatial extent of 
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the thermally activated transitions: Local thermal fluctuation can perhaps 
create small defects, but these topological defects will have a scale approx- 
imately equal to the size of their core (since both are defined by the same 
correlation length £). Such ill-defined defects are unlikely to survive — after 
all, they represent a configurations which may be a shallow local minimum 
of the free energy, but which have a higher free energy than the uniform 
superfluid. Hence, they are not large enough to be really topologically stable 
- change of the field configuration in a finite region of space the order of £ 
suffices to return to the uniform "true vacuum" . Moreover, such local defects 
cannot easily lead to the creation (or destruction) of the large ones. These re- 
quirements are least convincing in the case of monopoles (which have spatial 
extent of order £ in every direction). However even monopoles would have 
to be created "in pairs" by thermal fluctuations. These pairs of monopoles 
of opposite charge will be separated by distance of order £. Therefore, they 
are not likely to separate and survive. 

We are led to the conclusion that the dominant process in creation of 
topological defects will have to do with the critical slowing down and the 
subsequent "freeze-out" of the fluctuations of the order parameter at the 
time t rather than with the thermal activation and Ginzburg temperature. 

4 Numerical Experiments 

Our numerical experiment [14J to test the mechanism described in the pre- 
vious section for estimating the density of topological defects simulates the 
dynamics of a real field \& in a ID system. This system contains all the essen- 
tial ingredients needed to investigate topological defect formation scenarios 
(including symmetry breaking, although without an infinite range order), 
and at the same time, because of its one-dimensional nature, it allows us to 
perform a reasonable and accurate numerical analysis. We assume that the 
system under consideration is in contact with a thermal reservoir and obeys 
the Langevin equation 

* + 77* - d x J> + d*V(y) = 6 . (13) 

The potential V is given by ([!]). (A rescaling of x and t can be introduced 
to eliminate the parameters A and a.) The noise term has correlation 
properties 

(0(x, t), &(x', 0) = 2r]96(x' - x)5(t' - t) . (14) 

In Eq. (|14D, 77 is the overall damping constant characterizing the amplitude of 
the noise through Eq. (|13|). 9 fixes the temperature of the reservoir. Typical 
values used of 9 were irithe range (0.1, 0.01). 

Since the symmetry in the system is discrete (\1/ — > —^f), our numerical 
experiment consisted of investigating the creation of kinks (domain walls in 
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ID) as a function of the quench timescale tq. Figure shows a sequence of 
"snapshots" of \I/ for a typical quench. One can clearly see from this figure 
the transition from expectation value = above the critical temperature 
(e < 0) to = i^/e after the symmetry has been broken (e > 0). 

We test the predictions of defect density in [p]] in the overdamped regime 

(rj = 1), that is, the regime for which the damping term t]^/ in Eq. (|13|) 

is bound to dominate over For comparison, we also show in Figure |3] 
a sequence of "snapshots" of \I/ for an underdamped case with r] = 1/64. 
Overdamped evolution will take place whenever rf > |e| (see e.g. The 
characteristic relaxation time is then given by r = rj/\e\ , and therefore, from 
Eq. (H]), we have that fi = 1 and r « rj. Furthermore, from Eq. (^), one 
obtains that the freeze-out relaxation time in this case is 

i=VvT^- (15) 

One can then use (|T5|) to estimate from (||) and ( |T0|) both, the freeze-out 
temperature parameter e and correlation length £. Their values are 



e = i/r Q = y/rj/TQ (16) 

and 

t = to(T Q / v y/\ (I?) 

respectively, where we have adopted v = 1/2 in accord with the Landau- 
Ginzburg theory. The value of v = 1/2 is also in agreement with the calcu- 
lated critical exponents of the correlation length inferred from our equilib- 
rium simulations (see Fig. |j). We obtained for the equilibrium correlation 
length a best fit yielding £ = 1.38 ± 0.06, v = 0.41 ± 0.03 above T c , and 
£ = 1-02 ± 0.04, v = 0.48 ± 0.02 below T c , close to the Landau- Ginzburg 
exponent of v — 1/2. 

Figure | encapsulates the main result of our paper. It shows the number 
of zeros (kinks) obtained in a sequence of quenches with various values of the 
quench timescale tq. For each value of tq, 15 quenches were performed in a 
computational domain of size 2,048 units with periodic boundary conditions. 
This size of computational domain is large enough to avoid boundary effects. 
The number of kinks produced after the quench was obtained by counting 
the number of zeros of Above and immediately below Tc, there is a 
significant number of zeros which have little to do with the kinks (see Fig. 
|2|). However, as the quench proceeds, the number of zeros quickly evolves 
towards an "asymptotic" value (see Fig. ^j. This change of the density of 
zeros and its eventual "stabilization" is associated with the obvious change 
of the character of \l/ and with the appearance of the clearly defined kinks. 
By then, the number of kinks is nearly constant in the runs with long tq, 
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and, correspondingly, their kink density is low. Even in the runs with the 
smallest tq, there is still a clear break between the post-quench rates of the 
disappearance of zeros and the long-time, relatively small rate at which the 
kinks annihilate. (Kinks will eventually annihilate, but on slow timescales 
set by the thermal diffusion of kinks over distances of the order of their 
separation.) 

Our simulation clearly indicates |14| a scaling relation for the number 
density of kinks. We find a kink density (no. of zeros/2,048) of n = (0.087 ± 
0.007) Tq°' 28±0 ' 02 , when the kinks are counted at approximately the same 
£/tq value. By contrast, when the counting is done at the same value of 
absolute time t, n has a somewhat shallower slope (n oc Tq ' 24 ), reflecting 
higher annihilation rate at high kink densities. Within the estimated errors, 
this result is in agreement with the theoretical prediction of an scaling of 

n » (ti/tq)*. 

An important point to consider is that there are of course no "true" phase 
transitions in ID, and that instead of vortices or strings we are counting 
kinks. The remarkable agreement of our numerical results |TJj with the 



theoretical predictions suggests that the crucial ingredient in the scenario 
of Ref. |IJ is not that much the existence of a "mathematically rigorous" 
phase transition. The essential conditions required for this scenario are: the 
presence of symmetry breaking and the scalings (El) and (|B|) for the correlation 
length and relaxation time. 



5 Summary and Discussion 

The central issue of the work presented in this paper was a numerical test 
of the mechanism for defect formation in the course of second order phase 
transitions involving a non-conserved order parameter. Both laboratory [10, 
HI 12 1 and computer experiments [14] have added strong support to the 



proposal that the initial density of the topological defects will be set by the 
correlation length at the freeze-out instant —t; that is, at the moment when 
the relaxation timescale of the order parameter will be comparable to the 
characteristic timescale of the quench. At that instant, as a consequence 
of critical slowing down, perturbations of the order parameter become so 
sluggish that they cease to keep up with the equilibrium scalings, so that in 
the time interval [—t, i\ the order parameter cannot adjust to the changes of 
thermodynamic parameters induced by the quench. By the time t > t below 
Tc, the dynamics will "restart"; however, it will be too late for the system 
to get rid of the defects, by then the defects will be "set in concrete," that 
is, stabilized by the topology (Although subsequent evolution involving their 
interactions will ensue, and, generally alter their density). 
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Strong evidence that thermal activation is not responsible for the forma- 
tion of topological defects is provided by the experimental results obtained 
in the superffuid He 4 experiment by the Lancaster group []IIJ and by the 
He 3 experiments in Helsinki and Grenoble [[□], |12| . Together, they demon- 
strate that the argument in Ref. ffl is unaffectedby the thermal activation 
process invoked in Kibble's original discussion ([||, see also |Uj). Ginzburg 
temperature — that is, the temperature at which vacuum can be thermally 
"flipped" over the potential barrier in regions of size £ — is far below the 
critical temperature in the superfluid helium. 

Our numerical experiments also provide a strong confirmation of the the- 
oretical predictions of the non-equilibrium scenario |]]. The scaling of the 
kink density with tq follows closely theoretical expectations. These results 
are also supported by the dependence of the number of kinks on the damping 
parameter 77 (compare Fig. ^| and Fig. |3[), as well as by the preliminary re- 
sults of the computer experiments involving complex order parameter and/or 
more than ID space that we are currently investigating. 
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Figure 1: Characteristic behavior of the relaxation time (a) and the correla- 
tion length (b) during a quench. Dashed lines in (a) represent the symmetry 
breaking timescale \e/e\ = \t\. Dashed line in (b) represents the freeze-out 
correlation length during the time interval [—t, t\ . 
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Figure 2: Snapshots of $ during kink formation with a quench timescale of 
tq = 64 and damping parameter rj — 1. The figures, from top to bottom, 
correspond to i = -80, 14, 33, 51 and 301, respectively (compare with fOp. 
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Figure 3: Snapshots of \l/ during kink formation with a quench timescale of 
tq = 64 and damping parameter 77 = 1/64. The figures, from top to bottom, 
correspond to t — -80, 14, 33, 51 and 301, respectively. As expected from 
Eq. ( [I7D , the number of kinks is approximately half of what is seen in Fig. 0, 
where rj = 1 
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Figure 4: Equilibrium correlation length with solid lines representing a fit 
to f = £, \e\- u . The best fitting yields £ = 1-38 ± 0.06, v = 0.41 ± 0.03 
(X 2 = 1.2) above 7b, and £ = 1.02 ± 0.04, v = 0.48 ± 0.02 { X 2 = 3.7) below 
Tc, close to the Landau-Ginzburg exponent of v — 1/2. 
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Figure 5: Number of defects as a function of quench timescale. The number 
of kinks is calculated at a time £/tq = 4. The straight line is the best fit to 
N = N TQ a with a = 0.28 ± 0.02 and iV = 178 ± 14 ( X 2 = 1.96). 
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Figure 6: Correlation length £ and average number of zeros as a function 
of time in units of Tq; from bottom to top for £ and top to bottom for the 
average number of zeros, tq = 4, 8, 2048, 4096. The number of kinks used 
in Fig. [| were obtained at t/rg ~ 4, except in the large, computationally 
expensive, tq cases where an extrapolated value was used. 
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